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A SYSTEMATIC METHOD OF CALCULATING REDUCED MATRIX 
ELEMENTS OF SINGLE-PARTICLE OPERATORS 
by William F. Ford and Richard C. Braley 
Lewis Research Center 

SUMMARY 

The necessary mathematical apparatus for the evaluation of reduced matrix ele- 
ments of single- particle operators is introduced and developed, and the matrix elements 
are expressed as sums of integrals over the inelastic transition densities p ^g j(r). A 
biorthogonal set of functions called Gauss- Laguerre functions is then defined, in terms 
of which the densities are expanded. The reduced matrix elements are reexpressed as 
sums involving the coefficients pLSJn the transition density expansion. A number of 
simple applications to static and dynamic nuclear properties is presented, and the more 
complicated case of nucleon scattering form factors is treated in detail. Various prac- 
tical considerations related to the method are discussed and analyzed. Appendixes are 

devoted to the effects of distributed nuclear charge, the Coulomb form factor, and nu- 

20 

merical examples for neon-20 ( Ne) intended to serve as a test case for programs 
based on this method. 


INTRODUCTION 

Recent developments in the theory of nuclear structure have made it possible to de- 
scribe collective vibrational or rotational nuclear modes in terms of the motion of the 
individual nucleons. This opens the way to more rigorous tests of the microscopic de- 
scription, such as the calculation of elastic and inelastic cross sections, static nuclear 
properties, and electromagnetic transition rates. 

A great many of these calculations involve the evaluation of reduced matrix elements 
of single-particle operators between initial and final nuclear states. For some of the 
more complicated nuclear models, this can be a formidable task; and, furthermore, the 
details vary from model to model. 



If during evaluation of the matrix element, however, integration over the coordi- 
nates of the single-particle operator is left until last, the preceding integrations may be 
performed without any explicit reference to the single-particle operator. And, in fact, 
the result depends only on the initial and final states of the nucleus and on the nature of 
the transition involved. This leads directly to the concept of inelastic transition den- 
sities Plsj^’ essentially give, for the transition i — f, the probability for ab- 

sorption of angular momentum J = L + S by the last nucleon as a function of its distance 
from the nuclear origin. 

The idea of transition densities is not a new one. Such quantities were employed, 

for instance, by Haybron and McManus (ref. 1) in describing excitation of quadrupole 

12 16 40 

and octupole levels of carbon-12 ( *C), oxygen-16 ( O), and calcium-40 ( Ca) using 

the particle-hole model, and by Gillet and Melkanoff (ref. 2) in describing inelastic 
electron scattering from the same nuclei. More recently, Glendenning (ref. 3) has de- 
veloped a coupled- channel treatment of inelastic proton scattering using form factors 
which may be obtained from inelastic transition densities, and has emphasized that the 
coefficients needed for their calculation should be provided by the nuclear structure 
theorist. He gives examples for some currently popular nuclear models, and Braley 
and Ford (ref. 4) show how these coefficients may be constructed from projected 
Hartree-Fock wave functions. (These coefficients are essentially equivalent to the nu- 
clear overlap coefficients Sj(if[ab) of eq. (10).) Despite these occasional indications 
of interest, however, the convenience afforded by use of transition densities has been 
largely unnoticed. 

In this report we attempt to provide a systematic description of the use of inelastic 
transition densities in calculating single- particle matrix elements. We find that these 
calculations are greatly simplified if the transition densities are expanded in terms of 
a biorthogonal set of functions which we call Gauss- Laguerre functions. Each density 
is then specified by a number of coefficients, and evaluation of any matrix element re- 
duces to a simple sum over these coefficients. In order to be complete, we present all 
pertinent formulas, although some of these may be found elsewhere. Numerical examples 
for Ne are also presented to serve as a test case for programs based on our results. 


TRANSITION DENSITIES 


Consider the matrix elements of an operator £2 between initial and final nuclear 
states, 


^if = (JfM f | JjMj) 


( 1 ) 


2 



We define a generalized reduced matrix element of by 


<f IIS2 "i> JM (2) 

with M held fixed during the summation. (When f 2 transforms like the M L compo- 
nent of a spherical tensor of rank J, equation (2) reduces to the standard definition of 
ref. 5.) Inverting equation (2) yields 

n if =52< J i M i’ JM l J f M f ) < f " nili) JM (3) 

JM 

Here JM | J^M^) is a Clebsch-Gordan coefficient and J is shorthand for y2J + 1. 

Now suppose that is a sum of single- particle operators, 

A 

n= 22 0 ^ ( 4 ) 

n=l 


A J_1_ 

with x n s (r n , x ) representing the radial and spin-angle coordinates of the n 1 nucleon. 
We can isolate the dependence of the initial nuclear wave function on the coordinates of 

f Vi 

the n u nucleon by making the expansion 

| J lMi ) = A' : ‘/ 2 (, JM, i a m a I JjM,) * (x n )< M (j r . . if n .i? n+r ■ • ? A ) (5) 

am JM a 

a 



— m ~ 

(x) = <P R (r)& L ^ (x) 
a a 


( 6 ) 




22 ^ m r^ m sl jm ) Y Zm/^l/2m £ 


m t m s 


(7) 


is a member of an appropriately chosen complete set of single- particle functions. The 
reduced matrix elements of may now be expressed as 


3 



( 8 ) 


=E S j(tf |ab)(blif2 na) 

ab 


JM 


where the generalized single-particle reduced matrix elements are given by 


<b!K2"a) JM =£- Yj (j a m a’ JM l j b m b>< bm bl J2(x )l am a> 


*2 

•* b m a m b 


(9) 


and the nuclear overlap coefficients Sj(if |ab) by 



with U(abcd|ef) = efW(abcd; ef). 

j 2 m 2 I J 1 M 1> < J 1 M 1> j3 m 



In the derivation, use is made of the relation (ref. 

3 |jm> 


( 10 ) 

5) 


- ^ ] ^2 m 2’ ^3 m 3 ^l m l’ ^2^2 1*^) 3 ^1^2^ (H) 

J 2 M 2 

and the orthogonality of the Clebsch-Gordan coefficients. 

In most of the cases considered herein f2(x) may be written 

S2(x) = E "LSJM^^LSJ^ ( 12 > 

LSJM 

using the spin- angle tensors defined by 

■^LSjW= E <™ L ,SM s |jM)Y fr)a SM < 13 > 

m l m s L S 


4 



(The operator 

particles c?q = 
and spin- angle 


OgM * s a rank-S tensor in the spin-space of the particle; for 
S 

1 and Oj = 2S, where S is the spin operator. ) This permits 
integrations to be performed separately: 


spin-1/2 
the radial 


(bllftlla) JM -^(yi^LSjlPa) ^bl w LSJM 1^ ( 14 ) 

LS 


(The spin-angle reduced matrix element in equation (14) is the standard one of ref. 5. ) 

Examining equations (8) and (14), we see that it is natural to define a set of i — f 
transition densities 


p LSJ (r) j(if |ab)(j^ ) ||^ j gj||j a ) (p-^(r)cp^(r) (15) 

ab 


for then we have simply 


LS 

(fnfi" 1 ) =J^ sj( r ^LSJ M ( r )r 2 dr 


(16) 

(17) 


EXPANSION IN GAUSS-LAGUERRE FUNCTIONS 

To facilitate calculations involving the transition densities P^gj( r ), we expand them 
in terms of the very useful biorthogonal set of functions suggested by Sawaguri and 
Tobocman (ref. 6). These functions, which will be referred to as Gauss- Laguerre func- 
tions, have the following properties: 


jr L (a,/3r) = (/3r) L e"“^ r)2 if^ 1/2 (^ 2 r 2 ) 





a,Pr) 


(18) 

(19) 
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( 20 ) 


£ ^nL ( “’ ^ dr = 6 m, 

X) ^nL^’ ^nL ( “’ ^ = 5 (I> ~ r,) (21) 

n=0 r 

Here .2^(z) is a generalized Laguerre polynomial (ref. 7) or order c and index n. 

The Bessel transform of £r nl ^(a,Pr) is 

£ j L (kr)^- nL (Q!,/3r)r 2 dr = /3 3 C L (o?) ^ - — ^ ~ a,rk ) ( 22 ) 

C L («) = 2^ a 3 / 2+L/2 (1 _ a) -L/2 (23) 

y = [2j3\/a(l - a)]" 1 (24) 

and in the limit k - 0 this becomes 

/ r L !r nL (a,^r)r 2 dr = /3 _L (1 - 0 ,)“ L " 3 / 2 ^— (25) 

In terms of Gauss- Laguerre functions, the transition density may be written 

PLSJ^ = y^PLSJn^nL^’ ^ ^ 

n=0 

the coefficients being given by 

Pi*sjn=/ PLSJ< r ^nL(“-' Jr > r2clr <27> 

This leads directly to our basic result, 

< f " n "i) L S JM = ^P LSJ nX a, LSJM (r) ^'nL (0! ’^ r)r dr (28) 

n=0 u 


6 



The significance of equation (28) is that the radial integral involving and &~ nL 

may be done analytically, in closed form, and thus for any £2 the evaluation of 
(fn£2ni) LSJM is reduced to a relatively simple sum over the coefficients pj*g Jn - 


ILLUSTRATIONS 


Suppose, for instance, that we wish to calculate some static nuclear properties. 
First, it is necessary to check the wave function normalization, that is, to evaluate 

(JjM-IJjM.) = (inlni) 00 (29) 


Here the operator £2 is 


£2 = 1 = 




x )b^ 


X ^LSJ^I A 

n=l A n^l tatm L a 


6 L0 6 S0 


(30) 


Hence, 




L0 SO 


(31) 


Using equation (25) with L = 0 in equation (28), we find that 

<i" 1 Mi ) 00 




(32) 


It will turn out that pj*g Jn is nearly always accompanied by the factor [a/(a - l)] 11 , 
and so we define 


n 


-if 

^LSJn 


a; \ n if 


a - 1 


PlSJn 


(33) 
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Hence, 


(J.M i |j i M.) =^5(1 - a) - 3/2 


w . . 

E ~n 
p 000n 

n=0 


(34) 


Since 


differentiation of equation (25) with respect to a may be used to evaluate the integral 
involving 


co 


LSJM 


/ s V^4 n 2 k 

(r) = r 6 T a 5 


LO SO 


(35) 


From this we deduce the following expression for the mean square matter radius: 


‘ J i M i 




n=l 


= (l - a)" 5/2 

11/ 9 
Ap 


n 3\~ii 

H ~ 


a 2 


:J P OOOn 


(36) 


To obtain the electric quadrupole moment Q of the nucleus, we must evaluate 
(illf2 Ili) 2 o> where 


£2 = e 


X r n Y 20M 


'1 + T 


(37) 


n=l 


th 


and r is the z-component of the isospin operator for the n nucleon. Comparison 


n 


of equations (37) and (12) shows that for this case 

/16A 1/2 2/1 + t 


^LSJM^ " 6 


r I — — 1 6 L2 6 S0 


(38) 
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Since the operator (1 + t)/ 2 picks out the proton part PLgj( r ) °f the transition density 

P USJ^ = P 1SJ^ + P LSJ^ 

(the superscript if will be dropped for clarity when necessary), use of equation (25) 
with L = 2 in equation (28) yields 

Q.eK /2 ^ (1 -^£ ? P 02n (40, 

' 5 ' n=0 


To conclude this section, let us evaluate some simple dynamic nuclear properties. 
Here we are usually concerned with the absolute square of the matrix element n^, 
averaged over initial states and summed over final states: 


jf M.M f 


if 


J 2 JM J 2 


(fun ni) 


JM 


(41) 


Consider, for example, the reduced electromagnetic transition probability B(EX) 
for the i — f transition. Here the operator n is given by 


n = 



(42) 


and the sum over final states includes a sum (not shown above) over the photon orienta- 
tion jx. In this case, the reduced matrix element is independent of u, so that the sum 

2 2 

over ix gives a factor of X (=J ); use of equation (25) then yields 


B(EX) 


5 ? 


J3 -X (l 


a) -x- 3/2 £; 


o p 

p X0Xn 


n=0 


(43) 
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As an example of a case where Wlsj M (and hence the generalized reduced matrix 
element) depends explicitly on M, let us evaluate the Born approximation electron 
scattering form factor | F(q) | 2 , which is defined as above with 



(44) 


Here we have 


"lsjm^ 


— i L i T (qr)Y* 




6 so 6 lj 


(45) 


and so (using eq. (22)) 


F <9>| 2 = 42 ^ XX<)Ln f 

J 2 LM L 2 n=0 Z 


(46) 


The sum over M is readily carried out by using the addition theorem for spherical 
harmonics, yielding 


1 F(q) | 2 




^^LOLn^nl/ 1 " 


2 


(47) 


Equation (47) does not contain modifications which arise when the finite size of the pro- 
ton charge distribution is considered. These are treated in appendix A. 


NUCLEON SCATTERING FORM FACTOR 

A more complicated application occurs in the analysis of inelastic nucleon scatter- 
ing, where we need reduced matrix elements of (for example) a central nucleon- nucleon 
interaction 


10 


( 48 ) 


v*£ v «Vo> 

n=l 


where Xg denotes the projectile coordinates and 

1 

S=0 


V^n’^ =Z V s(l ? 0 ' ? nl) 7 sW * W 


(49) 


Following Tobocman, we introduce the Fourier transform of Vg: 

ik- r -ik- r 


/, A q /* ik- ? -ik- r_ 

V s(l F 0 “ F nl) = < 2jr ) J v S^ e e 


dk 


(50) 


If each plane wave is expanded in spherical harmonics, the angular integration over k 
is trivial and leads to 


v<V*-o> -ft E f v s( k )iL< kr o>iL( kr n» t2 * 

Wlsjm •'O 


(51) 


By comparing this with equation (12) we can identify ^gLjjyp l^en from equation (28) we 
have 


(fHV" 1 ) LgJM - (-1) L+S+J+M ^ I ^j( x 0 ) f LSj( r o) 


(52) 


where the nucleon scattering form factor f^gj^g) is given by 

00 

f LSj( r o) = (~ j) ^ p LS.TnJ^ r2 dr ^nl/“’ ^ ( 53 ) 

To evaluate f^gj(rg), we first carry out the integration over r by means of equa- 
tion (22): 


QO 

W r 0> =(-V 3c L (Q! ?S^LSJnj( ^il^O^nl/ 1 " «^k)k 2 dk 

W n=0 


(54) 
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Next suppose that Vg(r) can be expanded in Gauss- Laguerre functions, 


V S< r > = £ V Sm3W“'-< 3 ’ r > (55) 

m=0 

V Sm =jC”v s < r)^ n0 (a , ,(! ,r ) r2 dr (56) 

The choice L = 0 is appropriate because of the central nature of Vg. Again using 
equation (22), we find 


v Q (k) = 


4tt/3' 3 c 


0«”>E V Sm^o(‘ * a '’’ / ' k > 

m=0 


(57) 


This leaves us with 


-if — ~ 


f LSj( r o) “ C L(-) C Q(")^ ^LSJn V Sm 

n=0 m=0 


x f i L O jr 0 )^' nL (l " Q!,rk)^' m0 (l - a',y' k)k 2 dk (58) 

Jn 


The last step in the analysis is to observe that the product of the two Gauss- 
Laguerre functions may be written in the form 


^(i - or,*) jr m0 (l - a',y-k) = kV k ^ 

N=0 


(59) 


where 


— = 4a y 3 + 4a' , y'^ 
? 2 


(60) 
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i 



The coefficient is obtained by multiplying together the two series expansions for 

the Gauss- Laguerre functions, the result being 



t 

K=0 


The integration over k may now be carried out with the help of 



k 2N+L e -k /4? 2 j L(kro)k 2 ^ = 1 v ^ N: (2?) 2N+Lf3 ^ L (l,?r 0 ) 


(62) 


This yields an expression for f^jCrg) which, after some rearranging, may be written 


f LSj( r o) " £ f LSJN^Nl/ 0, 
N=0 


(63) 


where the form factor coefficients 


f LSJN are given by 



Although the series above look somewhat formidable, they are easily handled by the 
computer and, in practice, converge quite rapidly for short-range forces. (The ex- 
ceptional case of the long-range Coulomb force is treated in appendix B. ) Only the 

m = 0 term is necessary if V Q (r) has Gaussian shape, for instance, and only a few val- 
if ^ 

ues of PLgj n are nonzero if harmonic oscillator single- particle functions are used. 
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Even when the more realistic Wood-Saxon single- particle functions are used, it is 
demonstrated in appendix C that an accurate description of p^gj(r) is provided with 
about 20 terms. 


PRACTICAL CONSIDERATIONS 

From the foregoing analysis, there emerges a straightforward recipe for calculating 
matrix elements of single- particle operators from a knowledge of the nuclear wave func- 
tions: 

(1) The overlap coefficients Sj(if | ab) are obtained, and then, by means of equa- 
tion (15), the transition densities p^gj(r). 

(2) The densities are expanded in a series of Gauss- Laguerre functions, whose coef- 
ficients are obtained by means of equation (27). 

(3) The coefficients PLSJn may then be usec * in relatively simple formulas to ob- 
tain the desired matrix elements. 

In this section we shall discuss some of the practical aspects of this recipe. 

To begin with, the symmetry properties of the overlap coefficients can be simpli- 
fied by defining a related coefficient 


J’ 

Sj(if|ab) 

(65) 


for then 


Sj(fi|ba) = (-l) ja jb+Jf Ji Sj(if|ab) 


( 66 ) 


(The Sj are the nuclear structure coefficients discussed by Glendenning. ) 

The reduced matrix element of the tensor which is needed for constructing 

the transition densities by means of equation (15), is conveniently expressed in terms of 
the quantity 


„ab = \/47r 
S LSJ " 

J a 


i +(1/2) 

h) <yi^Lsjiiy 


(67) 
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which has the symmetry property 


f& - <-l> L * S4j gS J (68) 


When triangle inequalities are satisfied by the triads (Z & , L), (j a , j^, J), 

and the sum (£ & + Z b + L) is even, g^j is given by 

and (L, S, J), 

4oL = <la- J.ibjl J0 ) 

(69) 


(70) 

ab s LOL^ B t A ) 

L1Ltl jVa 

(71) 

where A = (L + J + l)/2, B = X a + X fc , and X = (l - j)(2Z + 1) - (1/2); otherwise g^j 
vanishes. 

Some care must be taken in generating the Gauss -Laguerre functions to avoid ex- 
cessive roundoff error. A simple, but highly accurate, scheme based on the recursion 
relation satisfied by the functions is as follows: Let = 1 and A Q = 0, and calculate 

S n = ( L - f) A n - (^n 

(72) 

^n+1 = ^n + A n + < n + ^n 

(73) 

A n + 1 = «Vl " ^n 

(74) 


Then 

2 

•P nL ( a >M = ((3r) L e~ a ^ T) cp n (75) 


To evaluate the sum of a series of Gauss- Laguerre functions, let <^ + j = = 0 

and calculate 


15 



(76) 


S n = A n + ( L + | A n + l-^n + l 

= ^n+1 + A n+1 + n ls n ^ 7) 

A n = ( Pn~ ^n+l < 78 > 

Then 

N 2 

2 A n *nL (a » /3r ) = (M L e~ a ^ S Q (79) 

n=0 


Similar formulas may be developed for the functions /3r). For large values of 

j3r, double-precision arithmetic may be required. 

It is desirable to calculate the coefficients pj^jn by direct numerical integration 
of p^gj(r) (although this could be avoided if harmonic oscillator single- particle functions 
were used), for the method is then independent of the choice of the set cpAr). The 

d. 

number of terms needed can be determined by comparing the original density P T ,p j( r ) 

with the truncated Gauss- Laguerre series representation. When harmonic oscillator 

if 2 

single-particle functions are used, of course, PLgj( r ) will be a polynomial in r times 
a factor of the form exp(-i'r^); hence, in this case, truncation error can be avoided 
completely by putting (1 - a)/32 = v. 

If nucleon scattering form factors f T ^ T (r n ) are desired, it is necessary to calculate 
the coefficients Vg m . When the potential has terms of Gaussian shape, the following 
integral may be used: 



r L e 


2 /p2 r 

•*mL< Qf, ’ /3,r)r ' 


dr 



(80) 


2 2-2 

with L = 0, where y = a'P' + R , Only the m = 0 term is nonzero if the par am - 

2 — 2 

eters are adjusted so that y = /3 T , that is, (1 - at')@' = R . If it is desired to employ 
several potentials with different ranges using only one parameter set, it is best to 
choose y = /3' for the largest range. 
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I 



II 


The Yukawa force 


V s (r) = V„ 


e -r/R 


(81) 


presents some problems, for good convergence of the series for Vg(r) cannot be ex- 
pected. Fortunately, all that is really required is that a' and /3 r be selected so that 
the series (eq. (57)) for Vg(k) is accurate enough for use in equation (54). The first two 
coefficients in this series are given by 


V S0 = 



- xe -x Erfc(x) 


(82) 



where x = (2/3'R and the rest satisfy 


(83) 



Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, April 13, 1970, 

129-02. 
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APPENDIX A 


EFFECTS OF DISTRIBUTED NUCLEON CHARGE 


When the operator £2(x) arises from the interaction between an electron and the 
nucleon at x, it should be replaced by 


ft c (x) = J' p c (r - rOfiCx'Jdr' 


(Al) 


where p^ represents the charge distribution of the nucleon. In this section, we shall 

indicate how this modification may be accomplished by simple changes in a,fi, and 
if 

^LSJn* 

Suppose, for example, that only the proton is to be considered charged, and its dis- 
tribution is to be Gaussian with an rms radius of R : 

r 


p c (r- r') = (2 Ti)' 






(A2) 


Then if S2(x) possesses the expansion (eq. (12)), it can easily be shown by expanding the 
plane waves in equation (A2) that 


LSJM 


(A3) 


where 



dk j L (kr)e 


~(l/6)k^Rp 



i L ( kr, ) a, LSJM^ r, ) r ' 2 dr ’ 


(A 4) 


(It is understood that w L S j M ( r ’) contains the factor (* + T )/ 2 -) 

Now in evaluating a reduced matrix element of by means of equation (28), we 
encounter the integral 


I = 



w LSJM (r) ^nL (0! ^ r)r2 dr 


(A 5) 
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p 

With the substitution of equation (A4) for the integration over r may be 

carried out by means of equation (22), leaving 



dk e 


"(l/6)k^Rp 


’^nL {1 “"’ rk) Jo 


,|2 


dr' 

(A6) 


Next we observe that 


-(l/6)k 2 R 2 

e P ^ L (1 " a > yk ) ~ a C ’ yk ) (A?) 


where 


R 


« c = « + 


6y* 


= “[ i+ fd'“)K) 2 ] 


(A 8) 


with y unchained. A second application of equation (22) (this time to the k- integration) 
then yields 


where 


1 = D nL f 0 “LSJM< r, »nL<“c- V )r ' 2 dr ' 


(A9) 




0 


,n+(L/2)/ 1 _ a yi+(L/2)+(3/2) 
D nL=|— ' ° 


V“c/ 


1 - a 


(A10) 


(All) 


Comparison of equations (A5) and (A9) shows that the effect of using in- 

stead of o)j^gj M (r) in evaluating I is equivalent to carrying out the latter evaluation 
with parameter values a ^ and /3^ and then inserting the factor D^. Consequently, 
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all subsequent formulas (such as eq. (46) for |F(q) | 2 , for example) remain valid if the 
replacements 


of — a. 


0 - 0 , 


if n if 
p LSJn nL^LSJn 


(A 12) 


are made throughout. 

It is instructive to apply this method to obtain the mean square nuclear charge ra- 
dius. Writing 



(A 13) 


we obtain (cf. eq. (36)) 


(J.MjofjjMj) 



(A 14) 


After the substitutions noted in equation (A 12) are made throughout, the resulting ex- 
pression simplifies to 


(J M |n c |j 1 M i > = (JjMjnlJjMj) + i^r 2 (l - «)' 5/2 |(^- i^?goon 

z ^ \ “ /n=0 


(A 15) 


To evaluate the sum over pg 00 n’ consider the following (cf. eq. (34)): 


1 = 




Z 


(1 - a)" 3 / 2 


y^PpQOn 

n=0 


(A 16) 
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Using this and equation (A 8) in equation (A 15) we find that 

(J i M.|S2 c |j.M i ) = + Rp 


(A 17) 


as expected. 
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APPENDIX B 


COULOMB FORM FACTORS 

For the long-range Coulomb force 



(Bl) 


the series derived earlier for nucleon scattering form factors converge only condition- 
ally, and a different method should be used. We begin with the familiar multipole ex- 
pansion of |r Q - r n | which leads to 

(fHVHi) jvgjM = Y LM^ r 0^ f L^ r 0^ 6 S0 6 LJ ( B2 ) 


f l/ r p) ~ : L ~^ ! e2 y^PLQLn r < r > L ^nL^ ,/3r ^ r dr 


where r^r^) is the smaller (larger) of (r Q , r). 
Consider first the integral 


W ' r 0 L_1 / ° r L ^ nL (»,?r)r 2 dr 


where 


<P n (x) = f Q e”^ 1_Q! ^ y »// B (y)dy 


(B3) 


(B4) 

(B5) 


(B6) 
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and 




(~ n;L + ? y ) 

r H) 


Now it can be shown (ref. 7) that 


/ 


= ^ 1 ( y ) = ^n^ y) ' ^n+l (y) 


Thus, integrating equation (B6) by parts gives 


<P n ( x ) = e ~^ “ )X »/ / n +1 ( x ) + ( x " “)[<P n ( x ) - <P n+1 ( x )] 


This two-term recursion relation for <p n (x) is readily solved, yielding 

n-1 






or - 1, 




Consequently, the integral I^rg) may be written 

n-1 


w = ferr 


,n 


1 ° <r ° ) ' S hr)- ^ 


The integral Ig(rg) cannot be put into elementary form, but a change of variable 
the limits of integration are 0 and 1 leads to 


V r (P = 


r 0 



(B7) 


(B8) 


(B9) 


(BIO) 


(BH) 
so that 


(B12) 
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n 

where x = (1 - Qf)(/3r 0 ) . For small x, the confluent hypergeometric function may be 
developed in a power series. For larger x, the following equation may be used: 


I 




Next we turn to the integral 


(B13) 



r ' L ' 1 3 ; - I 1 L <“-' !r > r2 


dr 




where now 

, n (x) =j f e ay ^(y)dy 

with 


^(y) = e- y ^ +(1/2) (y) 


(B14) 


(B15) 


(B16) 


(B17) 


In this case, 


^(y) 


satisfies 


/ 


^n( y ) dy = -^n" 1(y) = ^n-l (y) " ^n (y) 


and so proceding much as before we find that 




f a 
a - li 


vn 


® 0 < x ) 


k m 


ax . L- 1 , x 
e ^m ( x ) 


(B18) 


(B19) 
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In this case, <?q(x) can be easily evaluated, with the result 


K n< r 0> 



(B20) 
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APPENDIX C 


NUMERICAL EXAMPLES 

To illustrate parameter values for a typical case, we present some numerical re- 
20 

suits for Ne as examples of the various formulas given previously. These results are 
representative only, and are not intended as a serious effort to reproduce the experi- 
mental data. 

The nuclear wave functions were projected from axially symmetric intrinsic 

Hartree-Fock states with 20 orbits, employing for the most part harmonic oscillator 

single-particle functions spanning the Is, lp, 2s-ld, and 2p- If shells. From these wave 

functions the nuclear overlap coefficients Sj(if |ab) were determined by the methods of 

reference 4. These overlap coefficients, which are listed in table I, form the point of 

departure for the results given in the following paragraphs. 

Using a modified Simpson's rule algorithm with end-point corrections, with 100 in- 

if 

tervals of 0. 1 fm width, the density coefficients p T ^ Jr| were computed numerically for 

the elastic 0 + — 0 + transition and the inelastic 0 + — 2 + transition. The values 

a. = 0. 5593 and /3 = 0. 8 fm - * were chosen to match roughly the harmonic oscillator 

_o 

parameter v = 0. 282 fm . The density coefficients are listed in table II; nonzero 
values for n > 2 (LSJ = 000) and n > 3 (LSJ = 202, 212) are indicative of roundoff and 
truncation errors, but they are gratifyingly small. 

In table m we give the calculated and experimental nuclear matter radius and tran- 
sition rate B(E2), and in table IV the Born approximation electron scattering form 

o 

factors for selected values of q . For the nucleon scattering form factors, the 
Glendenning-Veneroni force (ref. 8) was used: 

2 

V(x) = -52. 0 e -(r/1 ’ 85) (P TE + 0. 6 P gE ) (Cl) 


Here we used 0' = 1 and adjusted a ' so that only V g Q was nonzero. 


The form factor 

coefficients are listed in table V. 

As a last example, we consider a Hartree-Fock wave function for ^ 8 Ne using Wood- 
Saxon single-particle functions spanning the Is, lp, and 2s-ld shells (ref. 9). No parti- 
cular attempt was made to search for optimum values of a and 0; for a typical choice, 
the transition density P^j( r ) and its truncated Gauss- Laguerre series representation 
shown in figure 1. (A larger number of terms is required than in the harmonic oscil- 
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lator case because the Wood-Saxon transition density falls off exponentially with dis- 
tance, whereas the individual terms in the series representation have a Gaussian en- 
velope. ) The correspondence is quite good, except for the regions where the densities 
are very small. 
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TABLE I. - OVERLAP COEFFICIENTS a FOR ELASTIC 


on 

AND INELASTIC TRANSITIONS IN “^Ne 


Elastic 


Inelastic transition, 

transition, b 


0 + - 

2 + 


0 + - 0 + 





a, b 

Sj(ii 1 ab) 

a, b 

Sj(ifjab) 

a, b 

Sj(if |ab) 

i,i 

1. 93009 

1,3 

-4. 48480 (-2) 

6, 10 

6. 74103 (-2) 

1, 2 

2.9913 (-1) 

1,4 

8. 53831 (-2) 

7, 6 

1. 40121 (-2) 

2,2 

7.01152 (-1) 

2,3 

1. 24356 (-1) 

7,8 

5. 23539 (-2) 

3, 3 

1. 29565 (-1) 

2,4 

-3. 92796 (-1) 

7,9 

8.93141 (-3) 

4, 4 

1. 23920 

3,1 

1.92589 (-2) 

8, 5 

-1. 16449 (-2) 

5, 5 

1. 95365 

3, 2 

-4. 50232 (-2) 

8, 6 

1. 56179 (-2) 

5, 7 

2. 39522 (-1) 

3, 3 

-3. 03809 (-2) 

8, 7 

-3.90190 (-3) 

6, 6 

3. 83722 

3,4 

-4. 93584 (-2) 

8,8 

7.07764 (-3) 

6, 8 

7. 49017 (-1) 

4, 1 

3. 83128 (-2) 

8,9 

-6. 94078 (-3) 

7,7 

3. 07065 (-2) 

4,2 

-1. 81644 (-1) 

8, 10 

1. 13321 (-2) 

8,8 

1. 49829 (-1) 

4, 3 

6. 20682 (-2) 

9, 5 

3. 34721 (-2) 

9, 9 

1. 83771 (-2) 

4, 4 

-3. 10634 (-1) 

9, 6 

2. 29125 (-2) 

10, 10 

1. 02197 (-2) 

5,6 

-3.03551 (-3) 

9,7 

3.77817 (-3) 



5,8 

2. 27451 (-2) 

9, 8 

4. 12532 (-3) 



5, 9 

7. 96833 (-2) 

9, 9 

-4. 24863 (-3) 



6, 5 

1.62547 (-3) 

9, 10 

-1. 02218 (-3) 



6, 6 

-5. 02013 (-3) 

10, 6 

3. 46070 (-2) 



6, 7 

-1. 37811 (-2) 

10, 8 

6. 02385 (-3) 



6, 8 

2. 18438 (-2) 

10, 9 

9. 47013 (-4) 



6, 9 

-3. 85836 (-2) 

10, 10 

-2. 59567 (-3) 


^he single-particle states <p are labeled by a = (1, . . . , 10) 
for (p a = ( ls 1 / 2 > 2 s i/ 2> ld 3/2’ ld 5/2’ lp l/2’ lp 3/2’ 

2p l/ 2 » 2 P3/2> lf 5/2’ lf 7/2^ 
b Sj(ii|ba) = Sj(ii [ ab) . 
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TABLE n. - DENSITY COEFFICIENTS 

FOR ELASTIC AND INELASTIC 
on 

TRANSITIONS IN Ne 

[Expansion parameters: a = 0. 5593; 
0 = 0. 8 fm" 1 .] 


Elastic Inelastic transition, 

transition, 0 + — 2 + 

0 + - 0 + 



Density coefficient 

n 

P000n 

P 202n 

^212n 

0 

5. 3423 (-1) 

5. 1295 (-2) 

-2.6324 (-3) 

1 

-2. 2250 (-1) 

-1. 2736 (-1) 

3.0126 (-3) 

2 

-5. 2752 (-4) 

-4.3476 (-3) 

1. 1925 (-4) 

3 

-4.6650 (-3) 

6.6777 (-7) 

-1. 8232 (-8) 

4 

6. 0398 (-7) 

-4.4295 (-9) 

2. 5105 (-10) 

5 

1.0555 (-8) 

5.4102 (-10) 

4. 5131 (-11) 

6 

4.7921 (-9) 

6.8837 (-9) 

-2. 3466 (-10) 

7 

1.0035 (-11) 

6.7184 (-9) 

-3. 2639 (-10) 

8 

-9. 5549 (-10) 

-1. 0776 (-9) 

-1.4097 (-10) 

9 

5. 2544 (-11) 

-9. 5644 (-9) 

1. 4020 (-10) 


TABLE III. - NUCLEAR MATTER RADIUS AND REDUCED 

ELECTROMAGNETIC TRANSITION 

on 

PROBABILITY FOR ^ u Ne 



Calculated 

Experimental 

Nuclear matter radius, 

2. 774 

2. 79 

r, fm 



Reduced electromagnetic transition 

160. 34 

286. 5 

probability, B(E2), e^fm^ 





TABLE IV. - BORN APPROXIMATION ELECTRON 

SCATTERING FORM FACTORS FOR ELASTIC 

90 

AND INELASTIC TRANSITIONS IN Ne 

Square of momentum 

Elastic transition, 

Inelastic transition, 

transfer, 

0 + - 0 + 

0 + - 2 + 

q 2 , 

fm' 2 

Form factor, | F(q) | 2 

0 

1.0 

0 

0. 5 

2.6856 (-1) 

6. 4541 (-3) 

1.0 

6. 5945 (-2) 

6. 8279 (-3) 

1. 5 

1. 3860 (-2) 

3. 5390 (-3) 

2. 0 

2. 1062 (-3) 

1. 1193 (-3) 

2. 5 

6.2199 (-4) 

1. 5903 (-4) 

3.0 

3. 1148 (-5) 

2. 0261 (-6) 

3. 5 

1. 3677 (-4) 

9. 0944 (-5) 

4. 0 

1.6683 (-4) 

1. 7703 (-4) 

4. 5 

1. 4272 (-4) 

2. 0258 (-4) 

5.0 

1.0329 (-4) 

1. 8223 (-4) 

5. 5 

6.7898 (-5) 

1. 4233 (-4) 

6.0 

4. 1961 (-5) 

1. 0117 (-4) 



TABLE V. - PROTON SCATTERING FORM 


FACTOR COEFFICIENTS f^^ FOR 
ELASTIC AND INELASTIC 
TRANSITIONS IN 20 Ne a 


[Expansion parameters: a = 0. 5593; 
(3 = 0. 8 fm' 1 .] 



Elastic transition, 
0 + - 0 + 

Inelastic transition, 
0 + - 2 + 


Form factor coefficients 

N 

f 000N 

f 202N 

f 212N 

0 

-3. 1382 (+3) 


-9. 0077 (+2) 

-4. 5644 

1 

1.0660 (+3) 


5. 8531 (+2) 

2. 2725 

2 

-8.6744 (+1) 


2. 5171 (+1) 

1.1510 (-1) 

3 

2.7019 (+1) 


1. 4187 (-2) 

4. 0569 (-5) 

4 

-3.0411 (-3) 


-2.6422 (-2) 

-8.4441 (-5) 

5 

-7. 7322 (-4) 


2.4913 (-2) 

7.7467 (-5) 

6 

4. 1770 (-4) 


-1. 5616 (-2) 

-4.6899 (-5) 

7 

-1. 2930 (-4) 


6.2597 (-3) 

1. 7943 (-5) 

8 

1. 8187 (-5) 


-1.4507 (-3) 

-3. 8944 (-6) 

9 

-7. 0469 (-7) 


1. 4755 (-4) 

3.6049 (-7) 


a See text for details of the nucleon-nucleon inter- 
action and its expansion. 
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■i m iii 


i hi min n ii 


(ziz - rsn 



Radius, r, fm 

Figure 1. - Effect of truncation on Gaijss-Laguerre series representation for a Wood 
Saxon inelastic transition density p[vj(r). Inelastic transition 0 + - 2 + for 
neon-20; expansion parameters: a=0.5; p=0.8fm _1 . 


NASA- Langley, 1970 24 


E-5366 



National Aeronautics and Space Administration 
Washington, D. C. 20546 

official business FIRST CLASS MAIL 



POSTAGE and fees pai 
national aeronautics 

SPACE ADMINISTRATION 


^ ^rs 70185 00903 

10U 001 ^caphniS LABORATORY /WL0L/ 

AIR FORCE WEAPONS m p X ! CO 8711.7 
.TRTLAMD AFB, Otw MfcXIGU 


ATT E. LOO BOWMAN, 


CHIEF, TECH. 


library 


POSTMASTHR: 


If Undeliverable ( Section 1- 
Postal Manual) Do Not Re 


''The aeronautical and space activities of the United Stales shall be 
conducted so as to contribute . . . to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.’’ 

— National Aeronautics and Space Act of 1958 

NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
ro merit NASA distribution in English. 

SPECIAL PUBLICATIONS : Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data* compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



